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Abstract 

The spin density matrix (SDM) used in atomic and molecular physics is 
revisited for nuclear physics, in the context of the radial density functional 
theory. The vector part of the SDM defines a "hedgehog" situation, which 
exists only if nuclear states contain some amount of parity violation. 

PACS: 21.10.-k, 21.10.Dr, 21.10.Hw, 21.60.-n 

1 Introduction 

The subject of density functionals (DFs) in nuclear physics [1, 2 is presently 
receiving intense attention. One of its difficulties is the handling of interactions 
that depend on spins. While there is a priori no theorem preventing a the- 
ory with simple densities from accommodating the influence of spin dependent 
forces, it is likely that a generalization of density profiles to "spin density" ones 
should, in practice, make the construction of a DF easier. The purpose of this 
paper is see whether one can adapt to nuclear physics the same concept [5J H] 
as that used for many years in atomic and molecular physics. 

Given the creation and annihilation operators, at^ and afrr, respectively, of 
a nucleon with spin ct = ±i at position r, and given a density operator D in 
many-body space, the SDM, is defined by its matrix elements, 

Paa' (r) = Tr ai^ af„, V . (1) 

In the following, we shall take advantage of the recent proof ^5;, based upon 
the rotational invariance of the nuclear Hamiltonian, that the nuclear DF is a 
scalar, namely a radial density functional (RDF); accordingly, it is understood 
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in the following, unless explicitly stated otherwise, that the density operator, 
T), in many-body space, is a scalar under rotations. Since practical calculations 
for a DF can eventually result in Kohn-Sham (KS) potentials [5], the approach 
described by the present paper, with its explicit treatment of spin, might give 
indications for the spin-orbit term in KS equations. 

The basic formalism for the SDM is explained in Sec. 2. A mandatory 
generalisation of the formalism is explained in Sec. 3. An illustrative example 
is provided in Sec. 4. We conclude in Sec. 5. 

2 Basic formalism 

We first relate the local creation and annihilation operators to those of an £s 
shell model, 

aL = ^"'^(^) ^/m(^) ailma ' "'P'^' = "^'i'^' (^) ^i' m' {f) ttn'e'm'a' ■ (2) 

Here the wave functions, fni{'>')Yim{f), represent the orbitals created by the 
operators a|j^^, with real radial form factors Lpni{r). The summation, X^n^m' 
runs over a complete basis of orbitals, assumed to be discrete for the sake of 
simplicity. A generalization with continuum orbitals brings no difficulty except 
for slightly less simple notations. Isospin labels are understood. 

We then rearrange the products, at^apcr', into their scalar and vector parts 
in spin space with the usual Clebsch-Gordan coefficients, 

Ar-SMs - E(-)'""'(^ ^ ^ ^ (^'\SMs) a^ar.', (3) 

aa' 

This gives, after inserting Eqs. ([2]), 

£m£'m' 

with 

Beml'm'SMsi^) = 

Y 'Pne{r)ipn'e'{r){-)^^'"'{^(T^ -a'\SMs}alf,^^^an'e'm'a' ■ (5) 

Next we recouple the orbital momenta carried by the operators Birntm' sms j 
Cu-LMSMsir) = {^TTii I' - m2\LM) B,^,i,^,sMs{r) , (6) 

SO that 

ApsMs = Y Vm{r)Ytm'{r){-f'-"'' E(^mf ~ m'\LM) CfjjLMSMs{r) . 

iml'm' LM 

(7) 
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Upon taking advantage of the relations i) between spherical harmonics, 

\0 Oj\m -to' fi J -'W'i' 
and ii) between Wigner 3j-coefRcients and Clebsch-Gordan ones, 



{Iml' ~m'\LM) ^{-y-^+^^ V2L+1 . 

'to —to 



£ £' 
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the orthogonality between Wigner 3j-coefficients, 
^(2L + 1) 

simplifies Eq. ^ into, 

E (-)^-*^rL-M(r)x 



LA 0-M f_i 



(8) 



(9) 



(10) 
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'LMSMsir)- (11) 



In Eq. PT|) . we used the facts that all numbers, ^, ^', L, Af, are integers and that 
ft i' L\ 

the 37'-coefficient, „ _ _ , vanishes unless I + I' + L is even. 
\0 0/ 

Finally, a recoupling of total orbital momentum and total spin yields, 
DLSJ^ir) - {LMSMs\Jti) x 

MMs 



EH" 



/(2^+l)(2^' + l) f££'L 







LMSMs 



(r) 



(12) 



so that. 
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In terms of D^sJtiif), the scalar or vector operators for the SDM now read 

A.sMs - E(-)^ ^LM{f) Y.{LMSMs\J,i) DLSJ^ir) ■ (14) 

LM Jfi 



With scalar density matrices T> in many-body space, there wih be vanishing 
traces, Tr Disj^(r) I?, unless J = /i = 0. In this case, the corresponding 
Clebsch-Gordan coefficient becomes, 

l\S+Ms 

{L M S Ms\0 0) - Sls 5m -Ms ^== , (15) 

V 2o + 1 

so that the SDM scalar or vector elements reduce to, 

{\Ms Y* (f) Y (r) 

Tr A,sMs -D = ^Jl'l Tr Dssoo{r) V = -^== Tr Dssm{r) V . 



(16) 



3 Generalization 



Two very different spin profiles emerge from the study made in Sec. 2. For the 
first of them, namely, for S* = 0, the result is simple, since, necessarily in this 
case, t and are equal, 

i^oooo(0 ^= Y\ Vne{r)iPn'e{r)al^^^^an'ema ■ (17) 
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For the second profile, i.e., for 5 = 1, spherical symmetry is ensured by the 
fact that all three spherical harmonics are multiplied by the same, radial form 
factor, which we denote Phh{f) in the following; we have a 'hedgehog" situation. 
Here we mean hedgehog-like in the sense that the vector spin field has only a 
radial dependency. It must be noticed, however, that only those pairs of particle 
orbital momenta {£,£'}, where \£ — £'\ < 1, can couple to L = 1. If ^ = £' , the 
f £ £' L\ 

3j-coefficient, ( ^ ^ ^ 1 , vanishes identically, since £ + £' + 1 becomes odd. 

Conversely, ii £ — £' = ±1, the corresponding products of operators, al^^^an'tm' , 
have an odd parity. Since parity violations in nuclear states are most often 
too tiny to be observable, the density operators V of interest always have an 
even parity. Therefore, if the traces, Tr C^^^i i _Ms i Ms do not vanish 

completely, then they will detect parity violations in P. A basic RDF, that uses 
Dqooo only, has no easy signature for parity violations. It is the occurrence of 
a tiny, but non- vanishing profile from Unoo that allows a more elaborate RDF 
theory to explicitly accommodate parity violations. 

For the sake of completeness, we show in Eq. (|18|) this "hedgehog" operator, 
Dhhir) = -Diioo/'v/S, the trace of which with V is the coefficient of liAfs(^) in 
Eq. p^ . It reads, upon taking advantage of Eqs. (O, (O, (IT^ and p^ . 

3v47r 

nn tniz m (T(T Ms 



£ £' 1 



V'(2^+l)(2£'-f 1) g {£m£' -m'\l -Ms) x 

^\^\^ ^^S) a^ma an't,n-<y' ■ (18) 
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A natural way to enlarge the theory to cases where the S = 1 form factor 
is not tiny consists in embedding the nucleus in an external field, Hi, that 
simultaneously breaks the rotational symmetry and the parity. To avoid loosing 
the advantage of an RDF, i.e., the reduction of three-dimensional calculations 
to one-dimensional ones, the symmetry breaking can be chosen as a minimal 
one, in the following way. Let Hi be a negative parity operator, bounded 
from below, that transforms as a vector under rotations. There is no need to 
assume that Hi is only made of local fields, Hi = J^i hi{fi,ai), where and ai 
denote the position and spin of the «th nucleon; any complicated Hi is allowed 
for the argument to come. What counts is that the extended Hamiltonian, 
H' = H + Hi, which is bounded from below, now contains, besides the basic 
scalar and positive parity H, a vector and negative parity component Hi . Then 
we use the "constrained search" definition [7] of a DF, 

F[p]=lniv^pTTH'V, (19) 

where now V is generalized into an arbitrary density operator, without symme- 
try properties. Here the symbol, V p, means that the minimization of the 
energy is performed over subsets in the 2? space that show a given spin density 
matrix p. Then the same argument, as that used in ^5 , to restrict 2? to be a 
rotation scalar, can be extended to restrict I? to be a mixture I?oi of a scalar 
and a vector. Next one can take advantage of Eq. ([Ti)) and derive F[p] from 
those few and radial profile operators, DLSJfj.i'''), where the conditions, S = 0,1 
and J = 0, 1 give limits to L via the usual triangular rules. 

To conclude this Sec, we note that a spin density DF is usually not very 
useful for an isolated nucleus, but becomes legitimate for a non-isolated one. 



4 Toy model for an illustrative example 

Consider a fictitious ^^O nucleus made of a full Os shell and an almost full Op 
shell and driven by a harmonic oscillator Hamiltonian, 

nlmar nljfir 

Here, temporarily, the isospin label, t = ±i, is explicit. The relation between 
is and jj creation operators (and, similarly, for annihilation ones) in this toy 
model reads, 

mcr j/i 

A Slater determinant, will describe this nucleus for our model. Assume that 
a perturbation of the harmonic oscillator slightly mixes the Opi orbitals with 
the Is 1 orbitals. The mixtures read, 

2 

I3\ 1 =cose6t, 1 -Hsineb'', i . (21) 

.7—2'/^''^ Op-,^,r ls^,/i,r ^ ' 



5 



We keep intact a core, made of the Os and Op| orbitals. The 16-body operator, 
T> = |^)(</>|, is still a scalar under rotations, but it has now a negative parity 
component at first order in e. Such a state, and similar density matrices, would 
justify the use of a "spin RDF" , with two profiles. 

Let |0) denote the fully closed Os and Op shells. At first order in e, the wave 
function under consideration is, = |0) + eX^r I''')' ^'^^^ 



\r) = {b^,,.,Jo,U. + ^isl-iJo,i.-lr) |0>. (22) 

Protons and neutrons will give equal matrix elements; hence, within an inessen- 
tial factor of 2, isospin labels and summations can again be omitted. Notice also, 
incidentally, that the particle-hole states |t), shown in Eq. (1^^ . do not repre- 
sent center-of-mass spurious shifts; the latter induce dipoles, not monopoles, in 
the one-particle-one-hole space. 

In the jj representation, we obtain for the S — case 

-Doooo(?-) = ^= 'Pni{r)(Pn'i{r) '^blejf.bn'ej^, (23) 



and the scalar profile, Tr£'oooo|'/')(0|7 has a vanishing contribution from the 
first-order matrix elements, (Ol-Dooook) = (t|-Doooo|0) = 0, because of the re- 
striction to equal values of £. The zeroth-order profile from the Os- and Op-shells, 
respectively, is obviously 



0000 1 



(24) 

with an inessential coefficient, ^2/7r, omitted for simplicity. 

Again for the jj representation, we find for the S — 1 (hedgehog) case. 

Tin' vJJ jj' rn'tn' a<j' Ms 

V(2^+l)(2f + 1) (^^ (fm/ -m'|l -Ms) x 

(i ^ i _ a'\\Ms) {lTa\a\3Vi) hi,^^ {£'m'^<j'\j' fi'} b„,,'r^' : (25) 
which reduces into, 

Dhh{r) - I] i'Y'^ V(2^ + 1)(2^' + 1) fndr) ^n'i'{r) 



47r . 

e e i\ i I I' 1 





where { } is a Wigner 6j symbol. The equalities, j — j' and ^ — jj! , refiect the 
fact that the LS coupling used in the previous section. Sec. 3, boils down to 
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Figure 1: Scalar profile of the toy model 
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Figure 2: Vector profile of the toy model 



total spin J = 0, as demanded by the scalar nature of the many-body density 
operator T). Accordingly, in a jj scheme, both the particle and the hole total 
spin labels must be equal. 

The zeroth-order matrix element in e that results from Eqs. ([^^ and 
(0|i?ft,/i|0), trivially vanishes. Upon a simple inspection of the first-order matrix 

bQpiiy |0), it is seen that the only non- vanishing 



elements, 



contributions come from the cases, {nijfj,} = Opiv and {n'i'jfi} = ^siv, be- 
cause of the restrictions on the values of £ and Here i' denotes the magnetic 
label of both the particle and the hole in Eq. ([25]). The two values of in |t) give 
the same contribution, similarly to the two isospin components. With a global 

omitted, the S = 1 form factor for 



factor, 4 ^S/tt 
the toy model reads 
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(27) 



Figures 1 and 2 show these scalar and vector profiles, respectively, in ar- 
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bitrary units to avoid unessential global coefficients and because we prefer to 
compare shapes. There is no need to stress how different their shapes are. 

5 Discussion 

We set out to investigate the possible role of the spin density matrix in the con- 
struction of the density functional for nuclei. Such spin densities have played an 
important role in atomic and molecular physics. However, the severe constraints 
of rotational invariance and parity for nuclei led to the result that the vector 
part of the spin density essentially vanishes in a nuclear DF that properly takes 
into accoimt such symmetries, namely, in an RDF. Thus, there is no way, in 
this approach, to explicitly describe spin properties in a nuclear RDF. On the 
other hand, the vector part becomes a signature of parity violation allowed in 
the RDF theory. We were able to legitimize the use of a spin density RDF, at 
the cost of introducing an external perturbation that has negative parity and 
transforms as a vector. Future studies are needed to understand the role of the 
spin-density-matrix formalism, when symmetries are broken by external forces. 
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